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Abstract 

Recent results for the coexistence of ferromagnetism and unconventional su- 
perconductivity with spin-tiplet Cooper pairing are reviewed on the basis of the 
quasi-phenomenological Ginzburg-Landau theory. New results are presented for 
the properties of phases and phase transitions in such ferromagnetic supercon- 
ductors. The superconductivity, in particular, the mixed phase of coexistence of 
ferromagnetism and unconventional superconductivity is triggered by the spon- 
taneous magnetization. The mixed phase is stable whereas the other supercon- 
ducting phases that usually exist in unconventional superconductors are either 
unstable, or, for particular values of the parameters of the theory, some of these 
phases are metastable at relatively low temperatures in a quite narrow domain 
of the phase diagram. The phase transitions from the normal phase to the phase 
of coexistence is of first order while the phase transition from the ferromagnetic 
phase to the coexistence phase can be either of first or second order depending 
on the concrete substance. The Cooper pair and crystal anisotropics are relevant 
to a more precise outline of the phase diagram shape and reduce the degener- 
ation of the ground states of the system but they do not drastically influence 
the phase stability domains and the thermodynamic properties of the respective 
phases. The results are discussed in view of application to metallic ferromagnets 
as UGe2, ZrZn2, URhGe. 

1. Introduction 

1.1. Notes about unconventional superconductivity 

The phenomenon of unconventional Cooper pairing of fermions, i.e. the formation 
of Cooper pairs with nonzero angular momentum was theoretically predicted pQ in 
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1959 as a mechanism of superfluidity of Fermi liquids. In 1972 the same phenomenon 
- unconventional superfluidity due to a p-wave (spin triplet) Cooper pairing of 3 He 
atoms, was experimentally discovered in the mK range of temperatures; for details 
and theoretical description, see Refs. [2113 E] Note that, in contrast to the standard 
s-wave pairing in usual (conventional) superconductors, where the electron pairs are 
formed by an attractive electron-electron interaction due to a virtual phonon exchange, 
the widely accepted mechanism of the Cooper pairing in superfluid 3 He is based on 
an attractive interaction between the fermions ( 3 He atoms) due to a virtual exchange 
of spin fluctuations. Certain spin fluctuation mechanisms of unconventional Cooper 
pairing of electrons have been assumed also for the discovered in 1979 heavy fermion 
superconductors (see, e.g., Refs. jHHEllZl) as well as for some classes of high-temperature 
superconductors (see, e.g., Refs. [HI El HOI HD IH 113 El 113 EH] ) • 

The possible superconducting phases in unconventional superconductors are described 
in the framework of the general Ginzburg-Landau (GL) effective free energy func- 
tional ^3] with the help of the symmetry groups theory. Thus a variety of possible 
superconducting orderings were predicted for different crystal structures [13 H3 HEl 
[23 EH E2| • A detailed thermodynamic analysis [HJHH] of the homogeneous (Meissner) 
phases and a renormalization group investigation [TTj of the superconducting phase 
transition up to the two-loop approximation have been also performed (for a three- 
loop renormalization group analysis, see Ref. [23; f° r effects of magnetic fluctuations 
and disorder, see [2UI23)- We shall essentially use these results in our present consid- 
eration. 

In 2000, experiments [2E] at low temperatures (T ~ 1 K) and high pressure (T ~ 1 
GPa) demonstrated the existence of spin triplet superconducting states in the metallic 
compound UGe 2 . This superconductivity is triggered by the spontaneous magnetiza- 
tion of the ferromagnetic phase which exists at much higher temperatures and coexists 
with the superconducting phase in the whole domain of existence of the latter be- 
low T ~ 1 K; see also experiments published in Refs. [211 [21] > and the discussion in 
Ref. |2H]- Moreover, the same phenomenon of existence of superconductivity at low 
temperatures and high pressure in the domain of the (T, P) phase diagram where the 
ferromagnetic order is present has been observed in other ferromagnetic metallic com- 
pounds (ZrZn 2 [HO] and URhGe [31] ) soon after the discovery j2EJ of superconductivity 
in UGe 2 . 

In contrast to other superconducting materials, for example, ternaty and Chevrel phase 
compounds, where the effects of magnetic order on superconductivity are also substan- 
tial (see, e.g., [221 E3 EH E3), in these ferromagnetic compounds the phase transition 
temperature (Tj) to the ferromagnetic state is much higher than the phase transition 
temperature {Tps) from ferromagnetic to a (mixed) state of coexistence of ferromag- 
netism and superconductivity. For example, in UGe2 we have Tps — 0.8 K whereas the 
critical temperature of the phase transition from paramagnetic to ferromagnetic state 
in the same material is Tf = 35K [23I2Z]- One may reliably assume that in such kind of 
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materials the material parameter T s defined as the (usual) critical temperature of the 
second order phase transition from normal to uniform (Meissner) supercondicting state 
in zero external magnetic field is quite lower than the phase transition temperature Tps- 
Note, that the mentioned experiments on the compounds UGe2, URhGe, and ZrZn2 do 
not give any evidence for the existence of a standard normal-to-superconducting phase 
transition in zero external magnetic field. 

Moreover, it seems that the superconductivity in the metallic compounds mentioned 
above, always coexists with the ferromagnetic order and is enhanced by the latter. As 
claimed in Ref. [26J in these systems the superconductivity seems to arise from the 
same electrons that create the band magnetism, and is most naturally understood as 
a triplet rather than spin-singlet pairing phenomenon. Note that all three metallic 
compounds, mentioned so far, are itinerant ferromagnets. Besides, the unconventional 
superconductivity has been suggested [36J as a possible outcome of recent experiments 
in Fe j37], in which a superconducting phase was discovered at temperatures below 
2 K at pressures between 15 and 30 GPa. Note, that both vortex and Meissner su- 
perconductivity phases [HZ] are found in the high-pressure crystal modification of Fe 
which has a hexagonal close-packed lattice. In this hexagonal lattice the strong ferro- 
magnetism of the usual bcc iron crystal probably disappears j3E]- Thus one can hardly 
claim that there is a coexistence of ferromagnetism and superconductivity in Fe but 
the clear evidence for a superconductivity is also a remarkable achievement. 

1.2. Ferromagnetism versus superconductivity 

The important point in all discussions of the interplay of superconductivity and ferro- 
magnetism is that a small amount of magnetic impurities can destroy superconductivity 
in conventional (s-wave) superconductors by breaking up the (s-wave) electron pairs 
with opposite spins (paramagnetic impurity effect |38j). In this aspect the phenomeno- 
logical arguments j3H] and the conclusions on the basis of the microscopic theory of 
magnetic impurities in s-wave superconductors jSH] are in a complete agreement with 
each other; see, e.g., Refs. jH2 HUE EH EH] • In fact, a total suppression of conventional 
(s-wave) superconductivity should occur in the presence of an uniform spontaneous 
magnetization M, i.e. in a standard ferromagnetic phase ^H]- The physical reason for 
this suppression is the same as in the case of magnetic impurities, namely, the opposite 
electron spins in the s-wave Cooper pair turn over along the vector M in order to lower 
their Zeeman energy and, hence, the pairs break down. Therefore, the ferromagnetic 
order can hardly coexist with conventional superconducting states. In particular, this 
is the case of coexistence of uniform superconducting and ferromagnetic states when 
the superconducting order parameter ip(x) and the magnetization M do not depend 
on the spatial vector x. 

But yet a coexistence of s-wave superconductivity and ferromagnetism may appear 
in uncommon materials and under quite special circumstances. Furthermore, let us 
emphasize that the conditions for the coexistence of nonuniform ("vertex", "spiral", 
"spin-sinosoidal" or "helical") superconducting and ferromagnetic states are less re- 
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strictive than that for the coexistence of uniform superconducting and ferromagnetic 
orders. Coexistence of nonuniform phases has been discussed in details, in particu- 
lar, experiment and theory of ternary and Chevrel-phase compounds, where such a 
coexistence seems quite likely; for a comprehensive review, see, for example, Refs. 

[321 E31 EH EH SO] 

In fact, the only two superconducting systems for which the experimental data allow 
assumptions in a favor of a coexistence of superconductivity and ferromagnetism are 
the rare earth ternary boride compound ErRl^B^ and the Chervel phase compound 
HoMoeSg; for a more extended review, see Refs. [SHUH]- In these compounds the phase 
of coexistence most likely appears in a very narrow temperature region just below the 
Curie temperature Tf of the ferromagnetic phase transition. At lower temperatures 
the magnetic moments of the rare earth 4/ electrons become better aligned, the mag- 
netization increases and the s-wave superconductivity pairs formed by the conduction 
electrons disintegrate. 

1.3. Unconventional superconductivity triggered by ferromagnetic order 

We shall not extend our consideration over all important aspects of the long stand- 
ing problem of coexistence of superconductivity and ferromagnetism rather we shall 
concentrate our attention on the description of the newly discovered coexistence of 
ferromagnetism and unconventional (spin-triplet) superconductivity in the itinerant 
ferromagnets UGe 2 , ZrZn 2 , and URhGe. Here we wish to emphasize that the main 
object of our discussion is the superconductivity of these compounds and, at a second 
place in the rate of importance we put the problem of coexistence. The reason is that 
the existence of superconductivity in such itinerant ferromagnets is a highly nontrivial 
phenomenon. As noted in Ref. [12] the superconductivity in these materials seems 
difficult to explain in terms of previous theories (221 EE1 EE] and seems to require new 
concepts to interpret the experimental data. 

We have already mentioned that in ternary compounds the ferromagtetism comes from 
the localized 4/ electrons whereas the s-wave Cooper pairs are formed by conduction 
electrons. In UGe2 and URhGe the 5/ electrons of U atoms form both superconductiv- 
ity and ferromagnetic order [2EJ EH] • In ZrZn 2 the same double role is played by the 4d 
electrons of Zr. Therefore the task is to describe this behavior of the band electrons at 
a microscopic level. One may speculate about a spin-fluctuation mediated unconven- 
tional Cooper pairing as is in case of 3 He and heavy fermion superconductors. These 
important issues have not yet a reliable answer and for this reason we shall confine our 
consideration to a phenomenological level. 

In fact, a number of reliable experimental data for example, the data about the coher- 
ence length and the superconducting gap [2S1 123 EH ED) , are in favor of the conclusion 
about a spin-triplet Cooper pairing in these metallic compounds, although the mecha- 
nism of this pairing remains unclear. We shall essentially use this reliable conclusion. 
Besides, this point of view is consistent with the experimental observation of coexistence 
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of superconductivity only in a low temperature part of the ferromagnetic domain of the 
phase diagram (T,P), which means that a pure (non ferromagnetic) superconducting 
phase has not been observed. This circumstance is also in favor of the assumption of 
a spin-triplet superconductivity. Our investigation leads to results which confirm this 
general picture. 

Besides, on the basis of the experimental data and conclusions presented for the first 
time in Refs. |26| and shortly afterwards confirmed in Refs. fI7\ I2E1 EDI EH] one 
may reliably accept the point of view that the the superconductivity in these magnetic 
compounds is considerably enhanced by the ferromagnetic order parameter M and, 
perhaps, it could not exist without this "mechanism of ferromagnetic trigger," or, 
in short, "iVf -trigger." Such a phenomenon is possible for spin-triplet Cooper pairs, 
where the electron spins point parallel to each other and their turn along the vector of 
the spontaneous magnetization M does not produce a break down of the spin-triplet 
Cooper pairs but rather stabilizes them and, perhaps, stimulates their creation. We 
shall describe this phenomenon at a phenomenological level. 

1.4. Phenomenological studies 

Recently, the phenomenological theory which explains the coexistence of ferromag- 
netism and unconventional spin-triplet superconductivity of Landau-Ginzburg type 
was developed [121 US]- The possible low-order couplings between the superconduct- 
ing and ferromagnetic order parameters were derived with the help of general symmetry 
group arguments and several important features of the superconducting vortex state 
in the ferromagnetic phase of unconventional ferromagnetic superconductors were es- 
tablished nam!. 

In this article we shall use the approach presented in Refs. |4*2} 14*3] to investigate the 
conditions for the occurrence of the Meissner phase and to demonstrate that the pres- 
ence of ferromagnetic order enhances the p-wave superconductivity. Besides, we shall 
establish the phase diagram corresponding to model ferromagnetic superconductors in 
a zero external magnetic field. We shall show that the phase transition to the su- 
perconducting state in ferromagnetic superconductors can be either of first or second 
order depending on the particular substance. We confirm the predictions made in 
Refs. |4*%1 H3] about the symmetry of the ordered phases. 

Our investigation is based on the mean-field approximation as well as on familiar 
results about the possible phases in nonmagnetic superconductors with triplet (p-wave) 
Cooper pairs [THl HH H2] • Results from Refs. (3111131311 will be reviewed and extended. 
In our preceding investigation |4^1 133 |3H] both Cooper pair anisotropy and crystal 
anisotropy have been neglected in order to clarify the main effect of the coupling be- 
tween the ferromagnetic and superconducting order parameters. The phenomenological 
GL free energy is quite complex and the inclusion of these anisotropics is related with 
lengthy formulae and a multivariant analysis which obscures the final results. 

Here we shall take into account essential anisotropy effects, in particular, the effect of 
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the Cooper pair anisotropy on the existence and stability of the mixed phase, namely 
the phase of coexistence of superconductivity and ferromagnetic order. We demonstrate 
that the anisotropy of the spin-triplet Cooper pairs modifies but does not drastically 
change the thermodynamic properties of this coexistence phase, in particular, in the 
most relevant temperature domain above the superconducting critical temperature T s . 
The same is valid for the crystal anisotropy, but we shall not present a thorough ther- 
modynamic analysis of this problem. The crystal anisotropy effect can be considered for 
concrete systems with various crystal structures [HI EE]- Here we find enough to demon- 
strate that the anisotropy is not crucial for the description of the coexistence phase. Of 
course, our investigation confirms the general concept ^H] that the anisotropy reduces 
the degree of degeneration of the ground state and, hence, stabilizes the ordering along 
the main crystal directions. 

There exists a formal similarity between the phase diagram obtained in our investiga- 
tion and the phase diagram of certain improper ferroelectrics ^TJ HHJ HHl GDI OHJ H2] ■ 
The variants of the theory of improper ferroelectrics, known before 1980, were criti- 
cized in Ref. for their oversimplification and inconsistency with the experimental 
results. But the further development of the theory has no such disadvantage (see, 
e.g., Ref. [HOI EI])- We use the advantage of the theory of improper ferroelectrics, 
where the concept of a "primary" order parameter triggered by a secondary order 
parameter (the electric polarization P e ) has been initially introduced and exploited 
(see Ref. j5UJ EH] E21)- The mechanism of the M-triggered superconductivity in itin- 
erant ferromagnets is formally identical to the mechanism of appearance of structural 
order triggered by the electric polarization P e in improper ferroelectrics (P-trigger). 
Recently, the effect of M-trigger has been used in a theoretical treatment of ferromag- 
netic Bose condensates 

1.5. Aims of the paper 

In the remainder of this paper we shall consider the GL free energy functional of uncon- 
ventional ferromagnetic superconductors. Our aim is to establish the uniform phases 
which are described by the GL free energy presented in Section 2.1. More information 
about the justification of this investigation is presented in Section 2.2. Note, as also 
mentioned in Section 2.2, that we investigate a quite general GL model in a situation 
of a lack of a concrete information about the values of the parameters of this model 
for concrete compounds (UGe2, URhGe, ZrZn 2 ) where the ferromagnetic superconduc- 
tivity has been discovered. On one side this lack of information makes impossible a 
detailed comparison of the theory to the available experimental data but on the other 
side our results are not bound to one or more concrete substances but can be applied to 
any unconventional ferromagnetic superconductor. In Section 3 we discuss the phases 
in nonmagnetic unconventional superconductors. In Section 4 the M-trigger effect will 
be described in the simple case of a single coupling (interaction) between the mag- 
netization M and the superconducting order parameter if) in an isotropic model of 
ferromagnetic superconductors, where the anisotropy effects are ignored. In Section 5 
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the effect of another important coupling between the magnetization and the supercon- 
ducting order parameter on the thermodynamics of the ferromagnetic superconductors 
is taken into account. In Section 6 the anisotropy effects are considered. In Section 7 
we summarize and discuss our findings. 

2. Ginzburg-Landau free energy 

Following Refs. [T%1 1^2*1 H3~] in this Chapter we discuss the phenomenological theory 
of spin-triplet ferromagnetic superconductors and justify our consideration in Sections 
3-6. 

2.1. Model 

Consider the GL free energy functional 

F[^M\ = J d 3 xf(^,M), (1) 

where the free energy density f(tp,M) (for short hereafter called "free energy") of a 
spin-triplet ferromagnetic superconductor is a sum of five terms: 

/(V, M) = / S (V0 + /f(M) + fity, M) + |J - B.M . (2) 

In Eq. (2) ip = = 1, 2, 3} is the three-dimensional complex vector describing the 

superconducting order and B = (H + AirM) = V x A is the magnetic induction; H 
is the external magnetic field, A = {Ay, j = 1,2,3} is the magnetic vector potential. 
The last two terms on the r.h.s. of Eq. (2) are related with the magnetic energy which 
includes both diamagnetic and paramagnetic effects in the superconductor (see, e.g., 

12211321152]). 

In Eq. (2), the term /s(^) describes the superconductivity for H = M = 0. This free 
energy part can be written in the form JH] 

/sW = hrM) + a s |Vf + |l# + yhAT + !E m 4 , (3) 

with 

fgradW = ^1 ( A^)*( A^) + ^2 [( A^)* + ( A^)* (^)] (4) 

+ir 3 (A^)*(A^), 

where a summation over the indices i,j (= 1, 2, 3) is assumed and the symbol 

d 2 1 e I 

A = -i/i7r- + —A' ( 5 ) 
ox, c 

of covariant differentiation is introduced. In Eq. (3), b s > and a s = a s (T — T s ), where 
a s is a positive material parameter and T s is the critical temperature of a standard 
second order phase transition which may take place at H = M. = 0; H = \H\, and 
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M. = \M\. The parameter u s describes the anisotropy of the spin-triplet Cooper pair 
whereas the crystal anisotropy is described by the parameter v s [TTJ[TE]. In Eq. (3) the 
parameters Kj, (j = 1,2,3) are related with the effective mass tensor of anisotropic 
Cooper pairs 



The term fp(M) in Eq. (2) is the following part of the free energy of a standard 
isotropic ferromagnet: 

&(M) =c f f2 |V,-M,| 2 + a/(T;)M 2 + |m 4 (6) 

3=1 

where Vj = d/dxj, bf > 0, and = ctf(T — Tf) is represented by the material pa- 

rameter atf > and the temperature Tf] the latter differs from the critical temperature 
Tf of the ferromagnet and this point will be discussed below. In fact, through Eq. (2) 
we have already added a negative term (— 2nM. 2 ) to the total free energy f(ip,M). 
This is obvious when we set H = in Eq. (2). Then we obtain the negative energy 
(— 2-kM 2 ) which should be added to f F (M). In this way one obtains the total free 
energy f-p(M) of the ferromagnet in a zero external magnetic field, which is given by 
a modification of Eq. (6) according to the rule 

f F (a f ) = f F [a f (T' f ) -> a f (T f )] , (7) 

where a,f = a/(T — Tf) and 

/ = '/ 



T f = r f + 2 -^ 



is the critical temperature of a standard ferromagnetic phase transition of second order. 
This scheme was used in studies of rare earth ternary compounds (221 EH ESI EE]- 
Alternatively [HZ] , one may work from the beginning with the total ferromagnetic free 
energy fp(af,M) as given by Eqs. (6) - (8) but in this case the magnetic energy 
included in the last two terms on the r.h.s. of Eq. (2) should be replaced with H 2 /87r. 
Both ways of work are equivalent. 

Finally, the term 

fi(^,M) =Z7 Q M.(V> x <0*)+ SM 2 \i/;\ 2 . (9) 

in Eq. (2) describes the interaction between the ferromagnetic order parameter M and 
the superconducting order parameter ip [121 H3|- The 7 -term is the most substantial 
for the description of experimentally found ferromagnetic superconductors |1H| while 
the <5ikf 2 1 ■?/> 1 2 -term makes the model more realistic in the strong coupling limit because 
it gives the opportunity to enlarge the phase diagram including both positive and 
negative values of the parameter a s . This allows for an extension of the domain of 
the stable ferromagnetic order up to zero temperatures for a wide range of values of 
the material parameters and the pressure P. Such a picture corresponds to the real 
situation in ferromagnetic compounds. 

In Eq. (9) the coupling constant 70 > can be represented in the form 70 = AnJ, 
where J > is the ferromagnetic exchange parameter [43J. In general, the parameter 
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S for ferromagnetic superconductors may take both positive and negative values. The 
values of the material parameters (T s , Tf, a s , a/, b s , u s , v s , bf, Kj, 70 and S) depend 
on the choice of the concrete substance and on intensive thermodynamic parameters, 
such as the temperature T and the pressure P. 

2.2. Way of treatment 

The total free energy (2) is a quite complex object of theoretical investigation. The 
various vortex and uniform phases described by this complex model cannot be investi- 
gated within a single calculation but rather one should focus on concrete problems. In 
Ref. the vortex phase was discussed with the help of the criterion [HBj for a stability 
of this state near the phase transition line T C 2(H); see also, Ref. [EH]- In case of H = 
one should apply the same criterion with respect to the magnetization M. for small 
values of \ip | near the phase transition line T c2 (A4) as performed in Ref. Here we 
shall be interested in the uniform phases, namely, when the order parameters ip and M 
do not depend on the spatial vector x G V (V is the volume of the superconductor). 
Thus our analysis will be restricted to the consideration of the coexistence of uniform 
(Meissner) phases and ferromagnetic order. We shall perform this investigation in de- 
tails and, in particular, we shall show that the main properties of the uniform phases 
can be given within an approximation in which the crystal anisotropy is neglected. 
Moreover, some of the main features of the uniform phases in unconventional ferro- 
magnetic superconductors can be reliably outlined when the Cooper pair anisotropy is 
neglected, too. 

The assumption of a uniform magnetization M is always reliable outside a quite close 
vicinity of the magnetic phase transition and under the condition that the supercon- 
ducting order parameter ip is also uniform, i.e. that vortex phases are not present at the 
respective temperature domain. This conditions are directly satisfied in type I super- 
conductors but in type II superconductors the temperature should be sufficiently low 
and the external magnetic field should be zero. Moreover, the mentioned conditions for 
type II superconductors may turn insufficient for the appearance of uniform supercon- 
ducting states in materials with quite high values of the spontaneous magnetization. 
In such cases the uniform (Meissner) superconductivity and, hence, the coexistence of 
this superconductivity with uniform ferromagnetic order may not appear even at zero 
temperature. Up to now type I unconventional ferromagnetic superconductors have not 
been yet found whereas the experimental data for the recently discovered compounds 
UGe2, URhGe, and ZrZn 2 are not enough to conclude definitely either about the lack 
or the existence of uniform superconducting states at low and ultra-low temperatures. 

In all cases, if real materials can be described by the general GL free energy (1) - 
(9), the ground state properties will be described by uniform states, which we shall 
investigate. The problem about the availability of such states in real materials at finite 
temperatures is quite subtle at the present stage of research when the experimental data 
are not enough. We shall assume that uniform phases may exist in some unconventional 
ferromagnetic superconductors. Moreover, we find convenient to emphasize that these 
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phases appear as solutions of the GL equations corresponding to the free energy (1) - 
(9). These arguments completely justify our study. 

In case of a strong easy axis type of magnetic anisotropy, as is in UGe 2 j2E], the 
overall complexity of mean-field analysis of the free energy f(ip, M) can be avoided 
by performing an "Ising-like" description: M = (0, 0, A4), where M. = ±|Af| is the 
magnetization along the "z-axis." Further, because of the equivalence of the "up" and 
"down" physical states (±Af) the thermodynamic analysis can be performed within 
the "gauge" M. > 0. But this stage of consideration can also be achieved without 
the help of crystal anisotropy arguments. When the magnetic order has a continuous 
symmetry one may take advantage of the symmetry of the total free energy f(ip, M) 
and avoid the consideration of equivalent thermodynamic states that occur as a result 
of the respective symmetry breaking at the phase transition point but have no effect on 
thermodynamics of the system. In the isotropic system one may again choose a gauge, 
in which the magnetization vector has the same direction as z-axis (\M\ = M z = A4) 
and this will not influence the generality of thermodynamic analysis. Here we shall 
prefer the alternative description within which the ferromagnetic state may appear 
through two equivalent "up" and "down" domains with magnetizations M. and (— M), 
respectively. 

We shall perform the mean-field analysis of the uniform phases and the possible phase 
transitions between such phases in a zero external magnetic field (H = 0), when the 
crystal anisotropy is neglected (v s = 0). The only exception will be the consideration 
in Sec. 3, where we briefly discuss the nonmagnetic superconductors (Ai ^ 0). For our 
aims we use notations in which the number of parameters is reduced. Introducing the 
parameter 

b = (b s + u s + v.) (10) 
we redefine the order parameters and the other parameters in the following way: 

(pj = b lf % = foe * , M = b) ,A M , (11) 
a s a f u s v s 

r = 7T i = 757' U ' = 1T' V = T' 

7o 8 

7 ~ b^bf ' 71 ~ (&M 1/2 ' 

Having in mind our approximation of uniform ip and M and the notations (10) - (11), 
the free energy density f(ip, M) = F(ip, M)/V can be written in the form 

f(ip,M) = r0 2 + i0 4 + 270 1 2 Msin(^-^ 1 )+ 7l 2 M 2 + tM 2 + iM 4 (12) 
-2w [0 2 2 sin 2 (# 2 - X ) + 2 0^sin 2 (0 1 - 6 3 ) + 2 2 sin 2 (# 2 - 6 3 )] 

Note, that in this free energy the order parameters if) and M are defined per unit 
volume. 
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The equilibrium phases are obtained from the equations of state 

df<JM>) 



dfi c 



0, (13) 



where the series of symbols fi can be defined as, for example, fi = {fi a } = (M, (pi, 3 , 
9i, 9 3 ); fi denotes an equilibrium phase. The stability matrix F of the phases /i is 
defined by 

F W = { F M] = . (14) 

An alternative treatment can be done in terms of the real (ipj) and imaginary (if)") 
parts of the complex numbers ipj = ipj + iip". The calculation with the moduli <f>j and 
phase angles 9j of ipj has a minor disadvantage in cases of strongly degenerate phases 
when some of the angles 9j remain unspecified. Then one should consistently use the 
properties of the respective broken continuous symmetry. Alternatively, one may do 
an alternative analysis with the help of the components ipj and if)". 

In order to avoid any ambiguity in our discussion let us note that we often use the term 
"existence" of a phase in order to indicate that it appears in experiments. This means 
that the phase, we consider, is either stable or metastable, in quite rare cases, when 
certain special special experimental conditions allow the observation of metastable 
states in equilibrium. When a solution (phase) of Eq. (13) is obtained it is said that 
the respective phase "exists", of course, under some "existence conditions" that are 
imposed on the parameters {fi a } of the theory. But this is just a registration of the 
fact that a concrete phase satisfies Eq. (13). 

The problem about the thermodynamic stability of the phases that are solutions of 
Eq. (13) is solved with the help of the matrix (14) and, if necessary, with an additional 
analysis including the comparison of the free energies of phases which correspond to 
minima of the free energy in one and the same domain of parameters {fi a }- Then 
the stable phase will be the phase that corresponds to a global minimum of the free 
energy. Therefore, when we discuss experimental situation in which some phase exists 
according to the experimental data, this means that it is a global minimum of the 
free energy, a fact determined by a comparison of free energies of the phases. If other 
minima of the free energy exist in a certain domain of parameters then these 

minima are metastable equilibria, i.e. metastable phases. If a solution of Eq. (13) is 
not a minimum, it corresponds to an (absolutely) unstable equilibrium and the matrix 
(14) corresponding to this unstable phase is negatively definite. 

When we determine the minima of the free energy by the requirement for a positive 
definiteness of the stability matrix (14), we are often faced with the problem of a 
"marginal" stability, i.e. the matrix is neither positively nor negatively definite. This 
is often a result of the degeneration of the states (phases) with broken continuous 
symmetry, and one should distinguish these cases. If the reason for the lack of a clear 
positive definiteness of the stability matrix is precisely the mentioned degeneration of 
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the ground state, one may reliably conclude that the respective phase is stable. If 
there is another reason, the analysis of the matrix (14) turns insufficient for our aims 
to determine the respective stability property These cases are quite rare and happen 
for very particular values of the parameters {fi a }- 

3. Pure superconductivity 

Let us set M = in Eq. (12) and briefly summarize the known results [TSl ITT] for the 
"pure superconducting case" when the magnetic order cannot appear and magnetic 
effects do not affect the stability of the normal and uniform (Meissner) superconducting 
phases. The possible phases can be classified by the structure of the complex vector 
order parameter ip = ("0i, "02, ^2)- We shall often use the moduli vector (0i, 2 , 03) 
with magnitude = (4>\ + 0| + ^l) 1 ^ 2 but we must not forget the values of the phase 
angles 9j. 

The normal phase (0,0,0) is always a solution of the Eqs. (13). It is stable for r > 0, and 
corresponds to a free energy / = 0. Under certain conditions, six ordered phases [T%]ITT| 
occur for r < 0. Here we shall not repeat the detailed description of these phases |18[ITT] 
but we shall briefly mention their structure. 

The simplest ordered phase is of type (-01, 0, 0) with equivalent domains: (0, ip 2 , 0) and 
(0, 0, Multi- domain phases of more complex structure also occur, but we shall not 
always enumerate the possible domains. For example, the "two-dimensional" phases 
can be fully represented by domains of type (-01, ip2, 0) but there are also other two 
types of domains: (-01, 0, ^3) and (0, -02, V^)- As we consider the general case when 
the crystal anisotropy is present {y 7^ 0), this type of phases possesses the property 

The two-dimensional phases are two and have different free energies. To clarify this 
point let us consider, for example, the phase (-01, ip2, 0). The two complex numbers, 
ipi and ip2 can be represented either as two-component real vectors, or, equivalently, 
as rotating vectors in the complex plane. One can easily show that Eq. (12) yields two 
phases: a collinear phase, when (82 — 61) = 7ik(k = 0, ±1,...), i.e. when the vectors 
ipi and ip2 are collinear, and another (noncollinear) phase when the same vectors are 
perpendicular to each other: (62 — 6\) = ir(k + 1/2). Having in mind that |0i| = 
|02| = 0/\/2, the domain (ipi, tp 2 , 0) of the collinear phase is given by (±1, 1, 0)</>/\/2, 
and the same domain for the noncollinear phase is given by (±i, 1, 0)0/ v^. Similar 
representations can be given for the other two domains of these phases. 

In addition to the mentioned three ordered phases, three other ordered phases exist. 
For these phases all three components ipj have nonzero equilibrium values. Two of 
them have equal to one another moduli <frj, i.e., <pi = <p2 = 4>3- The third phase is 
of the type 0i = 02 7^ 03 and is unstable so it cannot occur in real systems. The 
two three-dimensional phases with equal moduli of the order parameter components 
have different phase angles and, hence, different structure. The difference between 
any couple of angles 8j is given by ±7r/3 or ±27r/3. The characteristic vectors of 
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this phase can be of the form (e i7r / 3 , e~ i7r / 3 , l)<f>/y/3 and (e 2i7r / 3 , e" i27r / 3 , l)<f>/y/3. The 
second stable three dimensional phase is "real", i.e. the components if)j lie on the 
real axis; (9j — 9j) = 7fk for any couple of angles 9j and the characteristic vectors 
are (±1, ±1, 1)0/V3. The stability properties of these five stable ordered phases were 
presented in details in Refs. (HI UJ] • 

When the crystal anisotropy is not present (v = 0) the picture changes. The increase of 
the level of degeneracy of the ordered states leads to an instability of some phases and 
to a lack of some noncollinear phases. Both two- and three-dimensional real phases, 
where (9j — 9j) = 7rk, are no more constrained by the condition fa = <pj but rather have 
the freedom of a variation of the moduli <f>j under the condition 2 = — r > 0. The 
two-dimensional noncollinear phase exists but has a marginal stability [IT] . All other 
noncollinear phases even in the presence of a crystal anisotropy (v ^ 0) either vanish 
or are unstable; for details, see Ref. [TT] . This discussion demonstrates that the crystal 
anisotropy stabilizes the ordering along the main crystallographic directions, lowers the 
level of degeneracy of the ordered state related with the spontaneous breaking of the 
continuous symmetry and favors the appearance of noncollinear phases. 

The crystal field effects related to the unconventional superconducting order were es- 
tablished for the first time in Ref. [18J. In our consideration of unconventional ferro- 
magnetic superconductors in Sec. 4-7 we shall take advantage of these effects of the 
crystal anisotropy. In both cases v — and ti^O the matrix (14) indicates an instabil- 
ity of three-dimensional phases (all <pj ^ 0) with an arbitrary ratios 4>i/4>j- As already 
mentioned, for v ^ the phases of type fa = 02 7^ 03 are also unstable whereas for 
v = 0, even the phase fa = 02 = 03 > is unstable. 

4. Simple case of M-triggered superconductivity 

Here we consider the Walker- Samokhin model [13] when only the M0!0 2 — coupling 
between the order parameters ip and M is taken into account (7 > 0, 71 = 0). Besides, 
we shall neglect the anisotropies (w — v — 0). The uniform phases and the phase 
diagram in this case were investigated in Refs. ^J031I3S]- Here we summarize the main 
results in order to make a clear comparison with the new results presented in Sections 5 
and 6. In this Section we set 6*3 = and use the notation 9 = A9 = (0 2 — The 
symmetry of the system allows to introduce the notations without a loss of generality 
of the consideration. 

4.1. Phases 

The possible (stable, metastable and unstable) phases are given in Table 1 together 
with the respective existence and stability conditions. The stability conditions de- 
fine the domain of the phase diagram where the respective phase is either stable or 
metastable ^3]- The normal (disordered) phase, denoted in Table 1 by iV always ex- 
ists (for all temperatures T > 0) but is stable for t > 0, r > 0. The superconductivity 
phase denoted in Table 1 by SCI is unstable. The same is valid for the phase of co- 
existence of ferromagnetism and superconductivity denoted in Table 1 by C02. The 
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N-phase, the ferromagnetic phase (FM), the superconducting phases (SC1-3) and two 
of the phases of coexistence (COl-3) are generic phases because they appear also in 
the decoupled case (7 = 0). When the M0i0 2 -coupling is not present, the phases 
SC1-3 are identical and represented by the order parameter where the components 
(f>j participate on equal footing. The asterisk attached to the stability condition of 
"the second superconductivity phase" (SC2), indicates that our analysis is insufficient 
to determine whether this phase corresponds to a minimum of the free energy. As we 
shall see later the phase SC2, as well as the other two purely superconducting phases 
and the coexistence phase COl, have no chance to become stable for 7 7^ 0. This is so, 
because the non-generic phase of coexistence of superconductivity and ferromagnetism 
(FS in Table 1), which does not exist for 7 = is stable and has a lower free energy in 
their domain of stability. Note, that a second domain (M < 0) of the FS phase exists 
and is denoted in Table 1 by FS*. Here we shall describe only the first domain (FS). 
The domain FS* is considered in the same way. 

The cubic equation for M corresponding to FS (see Table 1) is shown in Fig. 1 for 7 = 
1.2 and t = —0.2. For any 7 > and t, the stable FS thermodynamic states are given 
by r(M) < r m = r(M m ) for M > M m > 0, where M m corresponds to the maximum 
of the function r(M). Functions M rn (t) and M (t) = (-t + 7 2 /2) 1/2 = V3M m (t) are 
drawn in Fig. 2 for 7 = 1.2. Functions r m (t) = 4M^(t)/7 for t < 7 2 /2 (the line of 
circles in Fig. 3) and r e (t) = 7|t| 1 ^ 2 for t < define the borderlines of stability and 
existence of FS. 



TABLE 1. Phases and their existence and stability properties [9 = (62 — #1), k = 0, ±1, ...]. 



Phase 


order parameter 


existence 


stability domain 


N 


tfij = M = 


always 


t > 0,r > 


FM 


4>j = 0, M 2 = -t 


t < 


r > 0, r 2 > 7 2 £ 


SCI 


01 = M = 0, 2 = — r 


r < 


unstable 


SC2 


02 = — r, 9 = irk, M = 


r < 


(t > 0)* 


SC3 


01 = 02 = M = 0, 4>l = -r 


r < 


r < 0, t > 


COl 


01 = 02 = 0, 2 = -r, M 2 = -t 


r < 0, t < 


r < 0, t < 


C02 


X = 0, 2 = -r, 9 = 9 2 = irk, M 2 = -t 


r < 0, t < 


unstable 


FS 


20f = 20 2 = 2 = -r + 7M, 03 = 


7M > r 


3M 2 > (-t + 7 2 /2) 




9 = 2vr(/c - 1/4), 7 r = ( 7 2 - 2t)M - 2M 3 




M > 


FS* 


20 2 = 20 2 = 2 = -(r + 7 M), 03 = 


—jM > r 


3M 2 > (-t + 7 2 /2) 




9 = 2n(k + 1/4), 7 r = (24 - 7 2 )M + 2M 3 




M < 



4.2. Phase diagram 

We have outlined the domain in the (t, r) plane where the FS phase exists and is a 
minimum of the free energy. For r < the cubic equation for M (see Table 1) and 
the existence and stability conditions are satisfied for any M > provided t > r ) 2 . For 
t < 7 2 the condition M > M have to be fulfilled, here the value M = (-t + 7 2 /2) 1/2 
of M is obtained from r(M ) = 0. Thus for r = the N-phase is stable for t > 7 2 /2, on 
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Figure 1: h — 77-/2 as a function of M for 7 = 1.2, and t = —0.2. 



M 



1- 
0.8- 




Figure 2: M versus t for 7 = 1.2: the dashed line represents M , the solid line represents 
M eq: and the dotted line corresponds to M m . 
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Figure 3: The phase diagram in the plane (t, r) with two tricritical points (A and B) 
and a triple point C; 7 = 1.2. The domains of existence and stability of the phases N, 
FM and FS are shown. 

the other hand FS is stable for t < 7 2 /2. For r > 0, the requirement for the stability 
of FS leads to the inequalities 

max M m ^j < M < M , (15) 

where M m = (M /a/3) and M should be the positive solution of the cubic equation 
of state from Table 1; M m > gives a maximum of the function r(M); see also Figs. 1 
and 2. 

The further analysis leads to the existence and stability domain of FS below the line AB 
given by circles (see Fig. 3). In Fig. 3 the curve of circles starts from the point A with 
coordinates (7 2 /2, 0) and touches two other (solid and dotted) curves at the point B 
with coordinates (— 7 2 /4, 7 2 /2). Line of circles represents the function r(M m ) = r m (t) 
where 



Dotted line is given by r e (t) = 7y \t\. The inequality r < r m (t) is a condition for the 
stability of FS, whereas the inequality r < r e (t) for (—t) > 7 2 /4 is a condition for 
the existence of FS as a solution of the respective equation of state. This existence 
condition for FS has been obtained from 7M > r (see Table 1). 

In the region on the left of the point B in Fig. 3, the FS phase satisfies the existence 
condition 7M > r only below the dotted line. In the domain confined between the 
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lines of circles and the dotted line on the left of the point B the stability condition for 
FS is satisfied but the existence condition is broken. The inequality r > r e (t) is the 
stability condition of FM for < {—€) < 7 2 /4. For (— t) > 7 2 /4 the FM phase is stable 
for all r > r e (t). 

In the region confined by the line of circles AB, the dotted line for < (— t) < 7 2 /4, 
and the t— axis, the phases N, FS and FM have an overlap of stability domains. The 
same is valid for FS, the SC phases and COl in the third quadrant of the plane (t, r). 
The comparison of the respective free energies for r < shows that the stable phase is 
FS whereas the other phases are metastable within their domains of stability. 

The part of the t-axis given by r = and t > 7 2 /2 is a phase transition line of second 
order which describes the N-FS transition. The same transition for < t < 7 2 /2 is 
represented by the solid line AC which is the equilibrium transition line of a first order 
phase transition. This equilibrium transition curve is given by the function 

r eq (t) = i [37 - ( 7 2 + let)) 1 / 2 ] M eq (t), (17) 

where 

M eq (t) = -L [ 7 2 - 8t + 7 ( 7 2 + 16t) 1/2 ] V2 (18) 

is the equilibrium value (jump) of the magnetization. The order of the N-FS transition 
changes at the tricritical point A. 

The domain above the solid line AC and below the line of circles for t > is the region 
of a possible overheating of FS. The domain of overcooling of the N-phase is confined 
by the solid line AC and the axes (t > 0, r > 0). At the triple point C with coordinates 
[0, r eq (0) = 7 2 /4] the phases N, FM, and FS coexist. For t < the straight line 



r*Jt) = 1 r + \t\, -7 2 /4<t<0, (19) 



describes the extension of the equilibrium phase transition line of the N-FS first order 
transition to negative values of t. For t < (— 7 2 /4) the equilibrium phase transition 
FM-FS is of second order and is given by the dotted line on the left of the point B 
(the second tricritical point in this phase diagram). Along the first order transition 
line r* (t) given by Eq. (fT9*|) the equilibrium value of M is M eq = 7/2, which implies 
an equilibrium order parameter jump at the FM-FS transition equal to (7/2 — 
On the dotted line of the second order FM-FS transition the equilibrium value of M is 
equal to that of the FM phase (M eq = y\t\). Note, that the FM phase does not exists 
below T s and this seems to be a disadvantage of the model (|T2*j) with ji = 0. 

The equilibrium phase transition lines of the FM-FS and N-FS transition lines in 
Fig. 3 can be expressed by the respective equilibrium phase transition temperatures 
T eq defined by the equations r e = r(T eq ), r eq = r(T eq ), r* q = r(T eq ), and with the 
help of the relation M eq = M{T eq ). This leads to some limitations on the possible 
variations of the parameters of the theory. For example, the critical temperature 
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{T eq = T c ) of the FM-FS transition of second order (7 2 /4 < —t) is obtained in the form 
T c = (T s + A<KJM/a s ), or, using M = {-a f /(3) l l\ 



where Tf > T s , and T* = (An J) 2 a f / a 2 (3 is a characteristic temperature of the model (fT^j) 
with 71 = w — v — 0. The investigation of the conditions for the validity of Eq. (|2*U|) 
leads to the conclusion that the FM-FS continuous phase transition (at 7 2 < — £) is 
possible only if the following condition is satisfied: 



where <j = (3a 2 /Aba 2 . This means that the second order FM-FS transition should 
disappear for a sufficiently large 7-coupling. Such a condition does not exist for the 
first order transitions FM-FS and N-FS. 

Taking into account the gradient term (4) in the free energy (j2J) should lead to a de- 
pression of the equilibrium transition temperature. As the magnetization increases 
with the decrease of the temperature, the vortex state should occur at temperatures 
which are lower than the equilibrium temperature T eq of the homogeneous (Meissner) 
state. For example, the critical temperature (T c ) corresponding to the inhomogeneous 
(vortex) phase of FS-type has been evaluated [HI] to be lower than the critical temper- 
ature (j2DI) : {T c — T c ) = AttubM-I 'o-s (a*b = \e\h/2mc - Bohr magneton). For J ^> fi B , 
we have T c m T c . 

For r > 0, namely, for temperatures T > T s the superconductivity is triggered by 
the magnetic order through the 7-coupling. The superconducting phase for T > T s 
is entirely in the (t, r) domain of the ferromagnetic phase. Therefore, the uniform 
supeconducting phase can occur for T > T s only through a coexistence with the ferro- 
magnetic order. 

In the next Sections we shall focus on the temperature range T > T s which seems 
to be of main practical interest. We shall not dwell on the superconductivity in the 
fourth quadrant (t > 0, r < 0) of the (t, r) diagram where pure superconductivity 
phases are possible in systems with T s > Tf (this is not the case for UGe2, URhGe, 
and ZrZn 2 ). Besides, we shall not discuss the possible metastable phases in the third 
quadrant (t < 0,r < 0) of the (t,r) diagram. 

4.3. Magnetic susceptibility 

Consider the longitudinal magnetic susceptibility xi — (Xv/^) per unit volume |4*fi] . 
The external magnetic field H = (0, 0, H) with H = (df /dM.) has the same direction 
as the magnetization M. We shall calculate the quantity x — yPfXi f° r the equilib- 
rium thermodynamic states fiQ given by Eq. (13). Having in mind the relations (11) 




(20) 



T f - T s > = ( ? + v/?)T* , 



(21) 
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between M and Ai, and between ip and y? we can write 



d 



X = 



dM n 



df_ 

dM 



(22) 



IH) 



where the equilibrium magnetization M and equilibrium superconducting order pa- 
rameter components ip j should be taken for the respective equilibrium phase (see 
Table 1, where the suffix "0" of 0, 9, and M has been omitted; hereafter the same 
suffix will be often omitted, too). Note that the value of the equilibrium magnetization 
M in FS is the maximal nonnegative root of the cubic equation in M given in Table 1. 

Using Eq. (22) we obtain the susceptibility x °f the FS phase in the form 



X' 1 = -l 2 + 2t + 6M 2 



The susceptibility of the other phases has the usual expression 

X' 1 = 2t + 6M 2 . 



(23) 



(24) 



Eq. (24) yields the known results for the paramagnetic susceptibility (xp = l/2t; 
t > 0) , corresponding to the normal phase, and for the ferromagnetic susceptibility 
(Xf — 1/4 j i j ; t < 0), corresponding to FM. These susceptibilities can be compared 
with the susceptibility x of FS. As the susceptibility x °f FS cannot be analytically 
calculated for the whole domain of stability of FS, we shall consider the close vicinity 
of the N-FS and FM-FS phase transition lines. 

Near the second order phase transition line on the left of the point B (t < — 7 2 /4), 
the magnetization has a smooth behaviour and the magnetic susceptibility does not 
exhibit any singularities (jump or divergence). For t > 7 2 /2, the magnetization is given 
by M = (s_ + s + ), where 



s± 



_jr ± 



(t ~ 7 2 /2) £ 
27 



- ' t)' 



1/2- 



1/3 



(25) 



For r = 0, M = 0, whereas for \^fr\ <C (t — 7 2 /2) and r = one may obtain M ps 
— 7r/(2t — 7 2 ) It. This means that in a close vicinity (r < 0) of r = along the 
second order phase transition line (r — 0,t > 7 2 ) the magnetic susceptibility is well 
described by the paramagnetic law xp — (l/2t). For r < and t — > 7 2 /2, we obtain 
M = -(7r/2) 1/3 which yields 



X 



= 6 



7 r 



2/3 



(26) 



On the phase transition line AC we have 

1 



M eq (t) = 



2V2 



7 



8t + 7 ( 7 2 + 16*) 



1/2' 



1/2 



(27) 
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and, hence, 

-v 2 r / 16A 1 / 2 ! 

28 



m\ 1/2 
3ll + 7j 



At the tricritical point A this result yields x -1 (v4) = 0, whereas at the triple point C 
with coordinates (0, 7 2 /4) we have x{C) = (2/7 2 ). On the line BC we obtain M = 7/2 
and, hence, 

X~ l = 2t + \ . (29) 
At the tricritical point B with coordinates (— 7 2 /4, 7 2 /2) this result yields x _1 (I?) = 0. 

In order to investigate the magnetic susceptibility tensor we shall slightly extend the 
framework of out treatment by considering arbitrary orientations of the vectors H and 
M. We shall denote the spatial directions (x,y,z) as (1,2,3). 

The components of the inverse magnetic susceptibility tensor 

xr 1 = rV&/ = te 1 } >A/ ( 3 °) 

can be represented in the form 

Xy 1 = 2(t + M 2 )5ij + 4MiMj + x ^\ , (31) 

where M and (pj are to be taken at their equilibrium values: M , <poj, 8 j. The last 
term in the r.h.s. of Eq. (28) is equal to zero for all phases in Table 1 except for FS 
(and FS*). When the last term in Eq. (29) is equal to zero we obtain the known result 
the susceptibility tensor for second order phase transitions (see, e.g., Ref. [T3j ). 

Consider the FS phase, where cj)j depends on Mj. Now we can choose again M = 
(0,0, M) and use our results for the equilibrium values of <f)j, 9 and M (see Table 1). 
Then the components Xij 1 corresponding to FS are given by 

X^ 1 = 2(t + M 2 )5ij + AMiMj - 7% . (32) 

Thus we have xi^j — 0, 

Xn=X 2 "2 1 = 2(t + M 2 ), (33) 
and X33 1 coincides with the inverse longitudinal susceptibility x" 1 given by Eq. (23). 
4.4. Entropy and specific heat 

The entropy S(T) = (S/V) = -Vd(f/dT) and the specific heat C(T) = (C/V) = 
T(dS/dT) per unit volume V are calculated in a standard way [TH]. We are interested 
in the jumps of these quantities on the N-FM, FM-FS, and N-FS transition lines. The 
behaviour of S(T) and C(T) near the N-FM phase transition and near the FM-FS 
phase transition line of second order on the left of the point B (Fig. 3) is known from 
the standard theory of critical phenomena (see, e.g., Ref. ^3] arid for this reason we 
focus our attention on the phase transitions of type FS-FM and FS-N for (t > — 7 2 /4), 
i.e., on the right of the point B in Fig. 3. 
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Using the equations for the order parameters ip and M (Table 1) and applying the 
standard procedure for the calculation of S, we obtain the general expression 

sm = -7=/ ~ > 2 ' (34) 

The next step is to calculate the entropies Sps{T) and Sfm of the ordered phases FS 
and FM. Note, that use the usual convention F N = V fx = for the free energy of the 
N-phase and, hence, we must set Sn(T) = 0. 

Consider the second order phase transition line (r = 0, t > 7 2 /2). Near this line 
Sps(T) is a smooth function of T and has no jump but the specific heat Cps has a 
jump at T = T s , i.e. for r = 0. This jump is given by 



ACW(T S 



1 



2t(T s )/ 7 2 J ' 



(35) 



The jump ACfs(T s ) is higher than the usual jump AC{T C ) = T c a 2 /b known from the 
Landau theory of standard second order phase transitions [T3] . 

The entropy jump ASac(T) = SpsiT) on the line AC is obtained in the form 



AS AC (T) = -M ( 



eq 



0£7 

4Vb 



16^ 1/2 



1+ , 

T 



^M eq } , (36) 



where M eq is given by Eq. (18). From Eqs. (18) and (36), we have AS(t = 7 2 /2) = 0, 
i.e., AS(T) becomes equal to zero at the tricritical point A. Besides we find from 
Eqs. (18) and (36) that at the triple point C the entropy jump is given by 

(37) 

On the line BC the entropy jump is defined by ASbc(T) = [Sps(T) — Sfm(T)]. We 
obtain 

A5ac(r) .( |t| _f)(- + ^). m 

At the tricritical point B this jump is equal to zero as it should be. The calculation of 
the specific heat jump on the first order phase transition lines AC and BC is redundant 
for two reasons. Firstly, the jump of the specific heat at a first order phase transition 
differs from the entropy by a factor of order of unity. Secondly, in caloric experiments 
where the relevant quantity is the latent heat Q = TAS(T), the specific heat jump 
can hardly be distinguished. 

4.5. Note about a simplified theory 

The consideration in this Section as well as in Sections 5 and 6 can be performed within 
an approximate scheme, known from the theory of improper ferroelectrics (see, e.g., 
Ref. |52j). The idea of the approximation is in the supposition that the order parameter 
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M is small enough so that one can neglect M 4 -term in the free energy. Within this 
approximation one easily obtains from the data for FS presented in Table 1 or by a 
direct calculation of the respective reduced free energy that the order parameters </> and 
M of FS are described by the simple equalities r = (7M — 2 ) and M = ( , y/2t)<f) 2 . Of 
course, one may perform this simple analysis from the very beginning. For ferroelectrics 
this approximation gives a substantial departure of theory from experiment [52]. In 
general, the domain of reliability of such an approximation should be the close vicinity 
of the ferromagnetic phase transition, i.e. temperatures near to the critical temperature 
Tf. On the other hand, this discussion is worthwhile only if the "primary" order 
parameter also exists in the same (narrow) temperature domain {(f) > 0). Therefore 
this approximation has some application in systems, where T s > Tf. 

For T s < Tf, one may simplify our thorough analysis by a supposition for a relatively 
small value of the modulus of the superconducting order parameter. This approx- 
imation should be valid in some narrow temperature domain near the line of second 
order phase transition from FM to FS. 

5. Effect of symmetry conserving coupling 

Here we consider the case when both coupling parameters 7 and 71 are different from 
zero. In this way we shall investigate the effect of the symmetry conserving 71-term in 
the free energy on the thermodynamics of the system. Note that when 7 is equal to 
zero the analysis is quite easy and the results are known from the theory of bicritical 
and tetracritical points ^3 HOI EM EH] ■ For the problem of coexistence of conventional 
superconductivity and ferromagnetic order this analysis (7 = 0, 71 7^ 0) was made in 
Ref. [321 • Once again we postpone the consideration of anisotropy effects by setting 
w = v = 0. The present analysis is much more difficult than that in Sec. 4, and cannot 
be performed only by analytical calculations; rather, some complementary numerical 
analysis is needed. Our investigation is based to a great extent on analytical calcu- 
lations but a numerical analysis has been also performed in order to obtain concrete 
conclusions. 

5.1. Phases 

The calculations show that for temperatures T > T s , i.e., for r > 0, we have three 
stable phases. Two of them are quite simple: the normal (N-) phase with existence 
and stability domains shown in Table 1, and the FM phase with the existence condition 
t < as shown in Table 1, and a stability domain defined by the inequalities r > jit 
and 



The third stable phase for r > is a more complex variant of the mixed phase FS and 
its domain FS*, discussed in Section 4. The symmetry of the FS phase coincides with 
that found in [131 

Let us also mention that for r < five pure superconducting (M = 0, <fi > 0) phases 
exist. Two of these phases, (0i > 0, (fi 2 = 4>3 — 0) an d (0i = 0, 2 > 0, 3 > 0) are 




(39) 



22 



unstable. Two other phases, (0i > 0, 2 > 0, 3 = 0, 9 2 = Q\ + nk) and (0i > 0, 2 > 
0,03 > 0,02 = Q\ + 7ck,8s - arbitrary; /c = 0, ±1,...) show a marginal stability for 
t > 7ir. 

Only one of the five pure superconducting phases, namely, the phase SC3, given in 
Table 1, is stable. In the present case of 71 7^ the values of <pj and the existence 
domain of SC3 are the same as shown in Table 1 for 71 = but the stability domain 
is different and is given by t > jir. When the anisotropy effects are taken into account 
the phases exhibiting marginal stability within the present treatment may receive a 
further stabilization. Besides, three other mixed phases (M 7^, <ft > 0) exist for r < 
but one of them is metastable (for jf > l,t < 71X, and r < 7it) and the other two 
are absolutely unstable. Here the thermodynamic behaviour for r < is much more 
abundant in phases than in the case of improper ferroelectrics with two component 
primary order parameter [SU]. However, at this stage of experimental needs about the 
properties of unconventional ferromagnetic superconductors the investigation of the 
phases for temperatures T < T s is not of primary interest and for this reason we shall 
focus on the relatively higher temperature domain r > 0. 

The FS phase is described by the following equations: 

X = 2 = A 3 = O , (40) 



2 = (± 7 M - r - 71M 2 ) , (41) 

(1 - 7l 2 )M 3 ± ^77iM 2 + - -L - 7ir ) M ± ^ = , (42) 

and 

(02 - 0i) = t\ + 2rrk , (43) 

(k = 0, ±1,...). The upper sign in Eqs. (41) - (43) corresponds to the FS domain in 
which sin(# 2 ~ &i) = —1 and the lower sign corresponds to the FS* domain. Here we 
have a generalization of the two-domain phase FS discussed in Section 4 and for this 
reason we use the same notations. The analysis of the stability matrix (14) for these 
phase domains shows that FS is stable for M > and FS* is stable for M < 0, just 
like our result in Section 4. As these domains belong to the same phase, namely, have 
the same free energy and are thermodynamically equivalent, we shall consider one of 
them, for example, FS. Besides, our analysis of Eqs. (40) - (43) shows that FS exists 
and is stable in a broad domain of the (t, r) diagram, including substantial regions 
corresponding to r > 0. 

5.2. Phase stability and phase diagram 

In order to outline the phase diagram (t, r) we shall use the information given above 
for the other three phases which have their own domains of stability in the (t, r) plane: 
N, FM, and FS. The phase diagram for concrete parameters of 7 and 71 is shown in 
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r 



1 




Figure 4: The phase diagram in the (t,r) plane for 7 = 1.2, 71 = 0.8 and w = 0. 

Fig. 4. The phase transition between the normal and FS phases is of first order and 
goes along the equilibrium line AC. It is given by the equation: 



eq 



(t) 



M, 



eq 



;i- 7l 2 )M e 2 g + ^77iM e(? +(t-^) 



(44) 



( 7 iM eg - 7/2) L v 

The equilibrium value M eq on the line AC is found by setting the equilibrium free energy 
/fs(/xo) °f FS equal to zero, i.e. equal to the free energy (f N = 0) of the N-phase. We 
have obtained the equilibrium energy f N as a function of the magnetization: 



M 2 



FS 



2(M 7l - 7/2) 2 
x /(l-7 1 2 )M 4 + 77 1 M 3 + 2 



(45) 



t{i-il)- 



r 



M 2 - 2 77l tM + t(t - — ) 



where M = M eq (hereafter the suffix "eq" will be often omitted). 

The numerical analysis of the free energy (45) as a polynomial of M shows that the 
expression in the curly brackets has one positive zero in the interval of values of t from 
t = 7 2 /2 (point A in Fig. 4) up to t = 0, where M t=0 = 7/2(71 + 1). As far as the 
obtained values for M are in the interval <M< (7/271) the existence condition of 
FS, namely, 

M(M 2 + 1) 



(7/2 - 71M) 



>0, 



(46) 



is also satisfied. 
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At the triple point C with coordinates t — 0, r = 7 2 /4(7! + 1) three phases (N, FM, 
and FS) coexist. To find the magnetization M on the equilibrium curve BC of the first 
order phase transition FM-FS for t< 0, we use the equality {fm=£fs, or, equivalently, 



(M 2 + t 2 ) 2 



l_M(l + 7l )] g + M(l- 7l ) 



2(M 7l - 7/2) 2 
Then the function r eq (t) for t < will have the form 



= 0. 



(47) 



4(1 + 71) 



(48 



This function describes the line BC of first order phase transition (see Fig. 4) which 
terminates at the tricritical point B with coordinates 



t B 



7 2 (2 + 7i) 



(49) 



4(1 + 7l ) 2 ' D 4(1 + 7l ) 2 ' 

To the left of the tricritical point the second order phase transition curve is given by 
the relation, 

r e (t) =lV^i + lit, (50) 

which coincides with the stability condition (39) of FM. This line intersects t-axis for 
t = (— 7 2 /7 2 ) and is well defined also for r < 0. On the curve r e (t) the magnetization 
is M = and the superconducting order parameter is equal to zero (0 = 0). The 
function r e (t) has a maximum at the point (t,r) = (— 7 2 /47 2 , 7 2 /47i); here M = 
(7/271). When this point is approached the second derivative of the free energy with 
respect to M tends to infinity, but as we shall see later the inclusion of the anisotropy 
of triplet pairing smears this singularity. The result for the curves r FS (t) of equilibrium 
phase transitions (N-FS ans FM-FS) can be used to define the respective equilibrium 
phase transition temperatures T FS . 

We shall not discuss the region, t > 0, r < 0, because we have supposed from the 
very beginning of our analysis that the transition temperature for the ferromagnetic 
order T/ is higher then the superconducting transition temperature T s , as i is for the 
known unconventional ferromagnetic superconductors. But this case may become of 
substantial interest when, as one may expect, materials with Tf < T s will be discovered. 

The stability conditions of FS can be written in the general form 



-M 2 + 77l M - t - 7 2 /2 
M71 - 7/2 



>0, 



(51) 



7M > , 



(52) 



1 



M71 - 7/2 



7l (l - 7i 2 )M 3 - | 7 (1 - 2 7l 2 )M 2 - ^ 7 2 7iM - l(t - 7 2 /2) 



> . (53) 
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Figure 5: The dependence M(t) as an illustration of stability analysis for 7 = 1.2 and 
71 = 0.8. 

Our consideration of the stability conditions (51) - (53) together with the existence 
condition Eq. (46) of the phase FS is illustrated by the picture shown in Fig. 5. 

For < t < 7 2 /2 and < M < (7/271) conditions (46) and (51) are satisfied. 
Condition (53) is a cubic equation in M(t) which for the above values of the parameter 
t has three real roots, one of them negative. The positive roots, M(t) > 0, as function 
of t are drawn by circles in Fig. 5 and it is obvious that the condition (53) will be 
satisfied for those values of M(t) that are between the two circled curves. The smaller 
positive root of Eq. (53) intersects t-axis for t = 7 2 /2 (point A in Fig. 5). Note, 
that M = 7/(271) is given by the horizontal dashed line. For t < — 7 2 (2 — 7 2 )/4 the 
stability condition (51) has two real roots shown by curves with crosses in Fig. 5. For 
negative values of the parameter t we shall consider also the curve M = \^—t which 
is the solution of existence condition (46) and is depicted by solid line in Fig. 5. For 
(-7 2 /47 2 ) < t < the FS phase exists and is stable when 7/(271) > M > y/^i. 

The point S in Fig. 5 with coordinates (— 7 2 / (47 2 ), 7/(271) is singular in sense that l.h.s. 
of conditions (51) and (53) go to infinity there. When t > (— 7 2 /47 2 ) the existence 
condition (46) implies 7/(271) < M < The stability condition (53) is always 

satisfied (two complex conjugate roots and one negative root) and condition (51) will 
be fulfilled for values of M between the two curves denoted by crosses in Fig. 5. 

5.3. Discussion 

The shape of the equilibrium phase transition lines corresponding to the phase tran- 
sitions N-SC, N-FS, and FM-FS is similar to that for the simpler case 71 = and we 
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shall not dwell on the variation of the size of the phase domains with the variations 
of the parameter 71 from zero to values constrained by the condition 7^ < 1 . Besides 
one may generalize our treatment (Section 4) of the magnetic susceptibility tensor and 
the thermal quantities in this more complex case and to demonstrate the dependence 
of these quantities on 71. We shall not dwell on these problems. But an important 
qualitative difference between the equilibrium phase transition lines shown in Figs. 1 
and 4 cannot be omitted. The second order phase transition line r e (t), shown by the 
dotted line on the left of point B in Fig. 1, tends to large positive values of r for large 
negative values of t and remains in the "second quadrant" (t < 0, r > 0) of the plane 
(t, r) while the respective second order phase transition line in Fig. 4 crosses the t-axis 
in the point t = — 7 2 /7i and is located in the third quadrant (t < 0,r < 0) for all 
possible values t < — 7V7 2 . This means that the ground state (at K) of systems 
with 71 = will be always the FS phase whereas two types of ground states, FM and 
FS, are allowed for systems with < 7 2 < 1. The latter seems more realistic in view 
of comparison of theory and experiment, especially, in ferromagnetic compounds like 
UGe2, URhGe, and ZrZn 2 . The neglecting of the 71-term does not allow to describe the 
experimentally observed presence of FM phase at quite low temperatures and relatively 
low pressure P. 

The final aim of the phase diagram investigation is the outline of the (T, P) diagram. 
Important conclusions about the shape of the (T, P) diagram can be made from the 
form of the (t, r) diagram without an additional information about the values of the 
relevant material parameters (a s , O/, •■■) and their dependence on the pressure P. One 
should know also the characteristic temperature T s , which has a lower value than the 
experimentally observed |2Hl EH EH1 EOl CHj phase transition temperature (Tps ~ 
to the mixed (FS) phase. A supposition about the dependence of the parameters a s 
and a/ on the pressure P was made in Ref . [IB] . Our results for Tf 3> T s show that the 
phase transition temperature Tp$ varies with the variation of the system parameters 
(ot a , ctf, ...) from values which are much higher than the charactestic temperature T s 
up to zero temperature. This is seen from Fig. 4. 

6. Anisotropy effects 

When the anisotropy of the Cooper pairs is taken in consideration, there will be not 
drastic changes in the shape the phase diagram for r > and the order of the respective 
phase transitions. Of course, there will be some changes in the size of the phase 
domains and the formulae for the thermodynamic quantities. The parameter w will also 
insert a slight change in the values of the thermodynamic quantities like the magnetic 
susceptibility and the entropy and specific heat jumps at the phase transition points. 

Besides, and this seems to be the main anisotropy effect, the w- and v -terms in the free 
energy lead to a stabilization of the order along the main crystal directions which, in 
other words, means that the degeneration of the possible ground states (FM, SC, and 
FS) is considerably reduced. This means also a smaller number of marginally stable 
states which are encountered by the analysis of the definiteness of the stability matrix 
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(14). All anisotropy effects can be verified by the investigation of the free energy (12) 
which includes the w- and f-terms. 



We have made the above general conclusions on the basis of a detailed analysis of the 
effect of the Cooper pair anisotropy (w-) term, as well as on the basis of a preliminary 
analysis of the total free energy (12), where the crystal anisotropy (v-) term is also 
taken into account. Here we shall present our basic results for the effect of the Cooper 
pair anisotropy on the FS phase; the crystal anisotropy is neglected (v — 0). 

The dimensionless anisotropy parameter w — u/{u + u) can be either positive or neg- 
ative depending on the sign of u. Obviously when u > 0, the parameter w will be 
positive too (0 < w< 1). We shall illustrate the influence of Cooper-pair anisotropy in 
this case. The order parameters (M, <pj, 9j) are given by Eqs. (40), (43), 

±7 M r - ^ 
(1-w) 

and 



3 

w - 7i)M 3 ± -77iM 2 + 
2 



7 2 

t(l — w) 7ir 



M ± — = , (55) 



where the meaning of the upper and lower sign is the same as explained just below 
Eq. (43). We consider the FS domain corresponding to the upper sign in the Eq. (54) 
and (55). The stability conditions for FS read, 

e-MhM-r-TjM* 
1 — w 

' 2W -("/M-wr- u»7iM 2 ) > 0, (57) 



and 



w 



1 — w 

I 2 

3(1 - w - 7i)M 2 + 3771M + t(l-w)--!-- -fir 



> . (58) 



For M 7^ (7/271 ) we can express the function r(M) defined by Eq. (54), substitute the 
obtained expression for r(M) in the existence and stability conditions (54)-(57) and do 
the analysis in the same way as for w — 0. The calculations show that in the domain 
r > 0, FS is stable for w < 0.5, when w = 0.5 there is a marginal stability, and for 
w > 0.5 the FS-phase is unstable (0 < w < 1). 

The results can be used to outline the phase diagram and calculate the thermodynamic 
quantities. This is performed in the way explained in the preceding Sections. 

7. Conclusion 

We have done an investigation of the M-trigger effect in unconventional ferromagnetic 
superconductors. This effect due to the Mipi ^-coupling term in the GL free energy 
consists of bringing into existence of superconductivity in a domain of the phase di- 
agram of the system that is entirely in the region of existence of the ferromagnetic 
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phase. This form of coexistence of unconventional superconductivity and ferromag- 
netic order is possible for temperatures above and below the critical temperature T s , 
which corresponds to the standard phase transition of second order from normal to 
Meissner phase - usual uniform superconductivity in a zero external magnetic field, 
which appears outside the domain of existence of ferromagnetic order. Our investiga- 
tion has been mainly intended to clarify the thermodynamic behaviour at temperatures 
T s < T < Tf, where the superconductivity cannot appear without the mechanism of 
M-triggering. We have described the possible ordered phases (FM and FS) in this most 
interesting temperature interval. 

The Cooper pair and crystal anisotropies have also been investigated and their main 
effects on the thermodynamics of the triggered phase of coexistence have been es- 
tablished. In discussions of concrete real material one should take into account the 
respective crystal symmetry but the variation of the essential thermodynamic prop- 
erties with the change of the type of this symmetry is not substantial when the low 
symmetry and low order (in both M and ip) 7-term is present in the free energy. 

Below the superconducting critical temperature T s a variety of pure superconducting 
and mixed phases of coexistence of superconductivity and ferromagnetism exists and 
the thermodynamic behavior at these relatively low temperatures is more complex 
than in known cases of improper ferroelectrics. The case Tf < T s also needs a special 
investigation. 

Our results are referred to the possible uniform superconducting and ferromagnetic 
states. Vortex and other nonuniform phases need a separate study. 

The relation of the present investigation to properties of real ferromagnetic compounds, 
such as UGe2, URhGe, and ZrZn 2 , has been discussed throughout the text. In these 
real compounds the ferromagnetic critical temperature is much larger than the super- 
conducting critical temperature (Tf ^> T s ) and that is why the M-triggering of the 
spin-triplet superconductivity is very strong. Moreover, the 71-term is important to 
stabilize the FM order up to the absolute zero (0 K), as is in the known spin-triplet fer- 
romagnetic superconductors. The neglecting [13] of the symmetry conserving 71-term 
prevents the description of the known real substances of this type. More experimental 
information about the values of the material parameters (a s , a/, ...) included in the free 
energy (12) is required in order to outline the thermodynamic behavior and the phase 
diagram in terms of thermodynamic parameters T and P. In particular, a reliable 
knowledge about the dependence of the parameters a s and aj on the pressure P, the 
value of the characteristic temperature T s and the ratio a s /a,f at zero temperature are 
of primary interest. 
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